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Abstract. 


Scattering  theory  is  developed  for  plane  diffraction  gratings. 

The  author's  theory  of  Rayleigh-Bloch  wave  expansions  is  used  to  construct 
wave  operators  and  a  scattering  operator  for  such  gratings.  For  gratings 
that  admit  no  surface  waves,  transient  wave  fields  near  the  gratings  are 
shown  to  behave  for  large  times  like  free  waves  and  corresponding 
asymptotic  wave  functions  are  calculated.  These  results  are  applied  to 
analyze  the  echoes  from  gratings  of  signals  due  to  localized  sources. 
Finally,  the  echoes  of  sources  remote  from  the  grating  are  estimated  and 
shown  to  be  completely  characterized  by  the  S-matrix  and  the  signal 
waveform. 


i 
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Introduction. 

A  theory  of  the  scattering  of  acoustic  and  electromagnetic  waves 
by  plane  diffraction  gratings  was  initiated  in  [5,  6].  That  work  developed 
the  theory  of  Rayleigh-Bloch  (or  R-B)  waves  for  diffraction  gratings. 
Physically,  the  R-B  waves  describe  the  response  of  a  diffraction  grating 
to  an  incident  monochromatic  plane  wave.  The  principal  result  of  [5,  6] 
is  an  R-B  wave  expansion  of  arbitrary  wave  fields.  The  purpose  of  the 
present  report  is  to  apply  the  results  of  [5,  6]  to  an  analysis  of  the 
scattering  of  transient  wave  fields  by  diffraction  gratings.  To  simplify 
the  analysis  it  is  assumed  throughout  this  report  that  the  grating 
propagator  A(G)  admits  no  R-B  surface  waves.  In  the  general  case  most 
of  the  results  of  this  report  hold  for  states  orthogonal  to  the  subspace 
spanned  by  the  surface  waves  (see  [5,  (6.36)]).  The  scattering  of  R-B 
surface  waves  is  not  analyzed  in  this  report. 

The  report  is  organized  as  follows.  §1  presents  a  review  of  the 
initial-boundary  value  problem  for  transient  wave  fields  near  a  diffrac¬ 
tion  grating  and  its  solution  by  means  of  the  Hilbert  space  theory  of 
the  grating  propagator  developed  in  [1,  3].  §2  presents  a  construction 

of  the  wave  operators  for  the  pair  consisting  of  the  reduced  propagator 
Ap  of  [5]  and  the  reduced  propagator  A0  p  for  the  degenerate  grating. 

The  results  of  §2  are  used  in  §3  to  construct  the  wave  operators  for  the 
pair  consisting  of  the  grating  propagator  A  and  the  propagator  Ag  of  the 
degenerate  grating.  The  results  of  §3  imply  that  transient  wave  fields 
near  gratings  have  asymptotic  wave  functions  in  the  sense  of  the  author's 
monograph  on  scattering  by  bounded  obstacles  [3].  These  functions,  and 
the  corresponding  asymptotic  energy  distributions  are  derived  in  54.  In 


.  * 


§5  Che  S-macrix  for  Che  pair  A  and  A„  is  calculaced.  16  concains  an 
analysis  of  Che  sCruccure  of  Che  echoes  produced  by  sources  chac  are  far 
from  che  gracing.  The  analysis  shows  chac  in  chis  case  Che  echo  wave¬ 
form  is  complecely  decermined  by  che  S-maCrix  and  che  signal  waveform, 
jusc  as  in  che  case  of  scaccering  by  a  bounded  obscade  [4]. 


Si.  The  Initial-Boundary  Value  Problem  for  the  Scattered  Fields. 

The  nocatlon  of  [5,  6]  will  be  used  in  this  report.  In  particular, 
X  *  (x,y)  €  R2  and  G  is  a  domain  in  R2.  Each  transient  acoustic  or 
electromagnetic  wave  field  in  G  can  be  described  by  a  real-valued 
potential  function  u  *  u(t,X)  that  is  a  solution  of  the  initial-boundary 
value  problem 

(1.1)  D2u  -  Au  «  0  for  all  t  >  0  and  X  €  G  , 

(1.2)  D^u  =  v  •  7u  *  0  (resp.,  u  -  0)  for  all  t  >  0,  X  €  3G  , 

(1.3)  u(0,X)  -  f(X)  and  Dtu(0,X)  -  g(X)  for  all  X  £  G  . 

Here  t  €  R  is  a  time-coordinate,  7u  •  (Du, Du),  Au  *  D2u  +  D2u,  3G 

x  y  x  y 

denotes  the  boundary  of  G  and  v(X)  is  a  unit  normal  vector  to  3G  at  X. 
u(t,X)  may  be  interpreted  as  a  potential  for  an  acoustic  field  with 
velocity  v  •  Vu  and  excess  pressure  p  *  Dtu.  With  this  interpretation  the 
boundary  condition  (1.2)  corresponds  to  an  acoustically  hard  (resp., 
soft)  boundary .  Alternatively,  if  u  satisfies  the  Neumann  condition  then 

(1.4)  E  -  Du  ,  E  -  -Du  ,  H  -  D.u 

x  y  y  x  z  c 

describes  a  TM  electromagnetic  field  in  a  domain  G  bounded  by  a  perfect 
electrical  conductor.  Similarly,  if  u  satisfies  the  Dirlchlet  boundary 
condition  then 

(1.5)  H  ■  -D  u  ,  H  *  D  u  ,  E  *  D  u 

x  y  y  x  z  c 


describes  a  TE  electromagnetic  field  in  the  same  kind  of  domain.  In 


c 
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both  che  acoustic  and  electromagnetic  interpretations  the  integral 


(1.6) 


E(u,K,t) 


{ (Dxu) 2 


(DyU)  1 


(Dtu) ‘ 


}  dxdy 


is  interpreted  as  the  wave  energy  in  the  set  K  at  time  t.  Note  that  all 
the  physically  observable  quantities  are  determined  by  the  derivatives 
of  the  potential  u.  The  state  of  the  wave  field  at  the  initial  time 
t  ■  0  is  defined  by  the  initial  conditions  (1.3). 

The  initial-boundary  value  problem  in  its  classical  formulation 
(1.1)  -  (1.3)  will  have  a  solution  only  if  3G  and  the  functions  f(X)  and 
g(X)  are  sufficiently  smooth.  However,  for  arbitrary  domains  G  the 
problem  is  known  to  have  a  unique  generalized  solution  with  locally 
finite  energy  whenever  the  initial  state  f,  g  has  this  property.  This 
result  was  proved  in  [1] .  In  cases  where  the  initial  state  has  finite 
energy. 


(1.7) 


f)2  +  (Dyf)2  +  g2)  dxdy  <  “  , 


a  simple  approach  to  the  initial -boundary  value  problem  is  provided  by 
the  grating  propagator 

(1.8)  A  -  AN(G)  (resp.,  A°(G)) 

of  [5,  §1].  For  arbitrary  domains  G  C  R2,  A  is  a  selfadjoint  realization 
in  the  Hilbert  space  L2(G)  of  the  operator  -A.  Moreover,  A  >  0  and 
D(A1,/2)  ■  Lj(G)  (resp.,  L2(G));  see  [5,  §1].  It  follows  that  if 


(1.9) 


f  e  L*(G)  (resp.,  L®(G))  and  g  €  L2(G) 
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chea  (1.7)  holds  and 

(1.10)  u(t,*)  -  (cos  t  Al/2)  f  +  (A"1'2  sin  t  A1/z)  g 

is  the  unique  solution  with  finite  energy  (*  solutions  wFE)  of  (1.1)  - 
(1.3);  see  [3,  5].  In  particular, 

(1.11)  u  €  C1 (R,L2 (G) )  n  C(R,D(A1/2  )) 

and  the  initial  conditions  hold  in  L2(G).  The  boundary  conditions  are 
incorporated  in  the  definition  of  D(A) .  The  d'Alembert  equation  (1.1) 
holds  in  a  suitable  weak  form;  see  [3].  The  transient  wave  fields 
studied  in  this  report  are  the  solutions  wFE  defined  by  (1.9),  (1.10). 

It  was  shown  in  [3,  Ch.  3]  that  solutions  wFE  in  arbitrary 
domains  have  a  representation 

1  /2 

(1.12)  u(t,X)  3  Re  {v(t,X) }  ,  v(t,*)  3  e  itA  h 

provided  that  f  and  g  satisfy  (1.9)  and  g  e  D(A-l/,z).  The  complex-valued 
function  h  €  D(A1',Z)  is  related  to  the  initial  state  f,  g  by 

(1.13)  h  *  f  +  i  A"I/Z  g  . 

This  representation  is  used  in  §4  below  to  determine  the  asymptotic 
behavior  for  t  08  of  the  transient  wave  fields  (1.10). 

The  R-B  wave  expansions  of  [5,  §6]  can  be  used  to  construct 
representations  of  the  solutions  wFE  (1.10)  and  (1.12).  In  the  case  of 
(1.12)  the  representations  take  the  form 
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(1.14)  v(t,X)  -  l.i.m.  f  (X,P)  e*itaj(P)  h+(P)  dP 

Jo2 

"•0 

where  P  ■  (p,q)  and  the  integral,  together  with  its  formal  t-derivative , 
converge  in  L2 (G)  [5,  Theorem  6.5  and  6.6]. 


§2.  Construction  of  the  Wave  Operators  for  A  and  A„ 
_ _ _ P  °,P 

The  notation  of  [5,  §3  and  §5]  is  used  in  this  section.  The 
purpose  of  the  section  is  to  prove  the  existence  and  completeness  of  the 
wave  operators 


(2.1) 


W  (A1^2  A1^2  J  ) 

vVAo,p’  P  ’  V 


itA?/2  -itAl/z 

s-lim  e  ,p  J,  -  p 
t-^±oo 


Q 


where  :  L2  (£2)  L2(Bfl)  is  defined  by 


(2.2) 


h(X) 


h(x)  ,  x  6  n  , 

0  ,  X6B,  -  SI  , 


This  will  be  done  by  means  of  an  explicit  construction  based  on  the 
eigenfunction  expansions  for  A^  and  A^  p  of  [5,  §3  and  §5].  The  principal 
results  are  formulated  as 

Theorem  2.1.  Let  G  be  a  grating  domain  [5,  §1].  Let 
p  €  (-1/2, 1/2]  and  assume  that  0„(Ap)  *  <J>.  Then  W+  p  and  W_  p  exist  and 
are  given  by 


(2.3) 


W 


*.P 


«L 

pP  +  fP 


In  particular,  W  :  L,  (fi)  -*■  L,  (B.)  are  unitary  operators  and  one  has 

—  »  P  W 

(2. A)  11^(1)  ■  w*  p  p(^)  w±  p  f°r  all  A  S  R  . 


Theorem  2.1  is  primarily  of  technical  interest  in  the  theory  of 


scattering  by  diffraction  gratings.  It  will  be  used  in  §3  to  derive  a 
construction  of  the  wave  operators  for  A  and  A„ . 
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Theorem  2.1  will  be  proved  by  the  method  of  [3,  Ch.  7].  Only 

the  case  of  W  will  be  discussed,  the  other  case  being  entirely 
■r»P 

similar.  To  begin  consider  the  wave  function 

-itA1/2 

(2.5)  v(t,‘)  =  e  p  h  ,  he  L2(0)  . 

The  eigenfunction  expansion  theorem  for  A^  of  [5,  Theorems  5.6  and  5.7] 
implies  that  v(t,X)  has  the  two  representations 

(2.6)  v(t,X)  *  l.i.m.  T  [  <£+  (X,p+m, q)  e“itw(P+m*9)  h+(p+m,q)  dq  , 

M^00  |m[ <M  •*  o  “ 

convergent  in  L2 (ft) .  As  in  [3,  Ch.  7],  the  incoming  representation  will 
be  used  to  calculate  the  behavior  of  v(t,*)  for  t  -► +°°.  The  eigenfunction 
<]>_  has  the  decomposition  [5,  (5.4)  and  (5.5)] 

(2.7)  <j>_(X,p+m,q)  *  j  (y)  (X,p+m,q)  +  <{>^(X,p+m,q) 

where  <K  is  incoming.  Combining  (2.6)  and  (2.7)  gives 

(2.8)  v(t,X)  *  j (y)  v*(t,X)  +  v+(t,X) 
where 

(2.9)  vJ(tfX)  -  l.i.m.  T  j  (X,p+m,q)  e-it^(P+m»9)  h_(p+m,q)  dq 

M-»o  |m|  <M  ‘  o 

converges  in  L2(B0)  while 

-M 

(2.10)  v  (t,X)  =  l.i.m.  7  $^_(X,p+m,q)  e  ^w(P*m’9)  (p+m,q)  dq 

M-*»  |m|<M  'o 

converges  in  L2 (ft) .  Note  that  the  convergence  of  (2.6)  and  (2.9)  implies 
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that  of  (2.10).  Moreover,  vQ(t,X)  is  a  wave  function  in  L2 (B  )  for  the 


reduced  propagator  A0  of  the  degenerate  grating;  namely. 


-itA 


1/2 


(2.11) 


vt(t,-)  =  e  ht 


where  h+  =  v+(0,‘)  e  L2(Bfl)  is  Siven  by 


(2.12) 


h0  =  $*  h  *  $  h 

0  o,p  -  o , p  -,p 


Theorem  2.1  will  be  shown  to  be  a  direct  corollary  of 

Theorem  2.2.  Under  the  hypotheses  of  Theorem  2.1  one  has,  for 
all  h  €  L2 (0) , 


(2.13) 

and  hence 


lim  v  (t,*)  *  0  in  L2  (ft) 
t-H-oo 


(2.14) 


lim  !v(t,*)  -  j(«)  v  (t , -) I  ,0.  =  0  . 
Proof  of  Theorem  2.2.  Equation  (2.14)  can  be  written 


(2.15) 


lim 

t-M-00 


-itA 


1/2 


t* 

-  J  e 


-itA 


1/2 

'’P  4> 

o,p  -,p 


hl  Lj  C£l)  •  0 


where  J  :  L2(B0)  -*■  L2(ft),  the  adjoint  of  the  operator  J  defined  in 

[5,  (5.22)],  is  given  by  J*  h(X)  *  j (y)  h(X)|^  .  Now  the  family  of 

operators  appearing  in  (2.15)  is  uniformly  bounded  for  all  t  €  R.  Hence 

to  prove  that  (2.15)  holds  for  all  h  S  L2 (ft)  it  will  suffice  to  verify 

it  for  all  h  from  a  dense  subset  of  L2 (0) .  It  will  be  convenient  to  use 

the  dense  subset  V~  ■  $*  V.  where 

o  -,p  o 


V0  C  l  9  L2(R0) 
m£Z 


(2.16) 
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is  the  set  of  all  g(q)  ■  tg^Cq)  :  m  €  z}  such  that  there  is  an  M  ■  M(g) 
with  the  properties 

(2.17)  gm(q)  -  0  for  lml  >  M  ,  and 

(2.18)  gm6  C>0  -  Em>p)  for  |m|  <  M 

where  is  the  exceptional  set  of  [5,  (5.15)].  Moreover,  it  will 

suffice  to  verify  (2.15)  for  functions  of  the  form 


(2.19) 


h(X) 


(X,p+m,q)  g(q)  dq 

0  ~ 


where  m  is  fixed  and  g  €  C.(Rn  -  E  )  has  support  in  an  interval 

u  u  m,  p 

I  C  Rfl  -  p,  since  the  case  of  a  general  h  S  Va  then  follows  by 
superposition.  Thus  the  proof  of  Theorem  2.2  may  be  completed  by  showing 
that  if 


(2.20) 


v+(t,X) 


4>'  (X,p+m,q)  e-ltuJ(P+m*9)  g(q)  dq 

J1 


where  g  €  C0(R  -  E  )  and  supp  g  c  I  then  (2.13)  holds. 
u  m ,  p 

The  definition  of  the  function  <}/ (X,p-Hn,q)  [5,  (5.4),  (5.5), 
(5.13),  (5.14)],  together  with  [5,  Theorem  4.15],  implies  that  for  fixed 
m  €  z  one  has 


(2.21)  <f^  6  C(fi  x  (r*  -  E)) 

where  E  is  the  exceptional  set  of  [5,  (2.30)].  Moreover,  the  far-field 
form  of  $1  is 
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L2 


(2.28)  llffl  §v^(t,*)lT  x  -  0  . 

C++®  L2(.Bo) 

It  remains  to  show  that  v£(t,*)  -*•  0  in  L2  (ft)  when  t  ■+■  +®.  This 
will  be  done  by  applying  (6,  Lemma  9.5]  to  u  ■  v2  *  v  -  v  .  To  this  end 
note  that  for  all  t  £  R  one  has  v+(t,*)  €  L2(ft)  by  (2.10)  and  v*(t,*) 

£  L2(B0)  by  [6,  (9.68)].  Thus  v*(t,*)  -  v+(t,*)  -  v*(t,«)  €  L2(ft)  for 
all  t  £  R  which  verifies  [6,  (9.75)]. 

The  local  decay  property  [6,  (9.76)]  follows  from  the  local 
compactness  property  of  the  grating  domain  G,  assumed  in  [5,  §1],  and 
the  abstract  decay  theorem  of  [2].  The  proof  is  the  same  as  that  of 
[3,  Theorem  5.5]  and  is  therefore  not  repeated  here. 

Finally,  [6,  (9.77)]  follows  directly  from  the  estimate  (2.25) 
and  the  representation 


(2.29) 


v+(t,X) 


P_(X,p+m,q)  e*ita,(P+m»<i)  g(q)  dq 


which  imply 


(2.30) 


|v* (t,X) 


<  C  e 


-uy 


|g(q) [  dq  for  all  X 


ft  , 

r 


and  t  €  R 


This  completes  the  proof  of  Theorem  2.2. 

Proof  of  Theorem  2.1.  The  proof  follows  that  of  [3,  Corollary 
7.2].  In  fact,  the  calculation  given  there,  adapted  to  the  present 
problem,  gives  the  estimate 
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(2.31) 


1  -It* ‘/! 

-itA1/2 

lh' 

0  » P  $ 

o,p  -,p 

h 

<  e 


P  -  J*  e 
* 


°’P  **  «  h 

«.P  ~,P  |L2(fl) 


-1CA  ' 

+  He  4,p<$*  4>  h)l.  .  +  ie 

0  » P  ~» P  1*2  (®n  y ' 


L,(«  ) 

2  o,r 


The  first  term  on  the  right  in  (2.31)  tends  to  zero  when  t  ■*  +»  by  (2.15) 
The  last  two  terms  tend  to  zero  when  t  +00  by  the  local  decay  property 
used  in  the  proof  of  Theorem  2.2.  It  follows  that  the  left-hand  side  of 

(2.31)  tends  to  zero  when  t  ♦  +00  which  proves  the  existence  of  W  and 

£  «P 

equation  (2.3).  Finally,  to  verify  (2.4)  note  that  it  can  be  written 

(2.32)  n  (X)  n  (X)  $*  *.  for  X  €  R 

P  i,P  o,P  o»P  o.p  ±,p 

by  (2.3).  The  unitarity  of  implies  that  an  equivalent  relation  is 

-  tP 


(2.33) 


„  Hn (X)  ♦*  -  $  IT  (X)  **  for  X  €  R  . 

p  ifp  o ,p  o ,p  0 ,p 


But  this  last  equation  is  correct  because  the  two  sides  coincide  with 
the  operation 


(2.34)  {gm(q)>  (H(X  -  U)2  (p+m,q) )  gm(q)} 

in  Z  ©  L2(R0);  see  [5,  Theorem  5.7].  This  completes  the  proof. 


t 


I 


f 


f 


» 


§3.  Construction  of  the  Wave  Operators  for  A  and  An. 

The  notation  of  [5,  §6]  Is  used  In  this  section.  The  purpose  of 
the  section  Is  to  prove  the  existence  and  completeness  of  the  wave 
operators 


(3.1) 


W±  -  W+(A01/2  ,Al/2  ,JQ) 


s-lim  e 


ItAj'1 


JG« 


-itA 


1/2 


where  :  L2  (G)  ■*  L2(Rj)  Is  defined  by 


(3.2) 


JG  h (X) 


h(X)  ,  X  6  G  , 

0  ,  X  e  Rj  -  G  . 


The  principal  results  of  the  section  are  formulated  as 

Theorem  3.1.  Let  G  be  a  grating  domain  [5,  §1]  and  let  A  »  A(G) 
admit  no  surface  waves.  Then  tf+  and  W_  exist  and  are  given  by 

(3.3)  -  ♦*  *_  • 

£  o 


In  particular,  W+  :  L2 (G)  L2(R2)  are  unitary  operators  and  one  has 

(3.4)  11(A)  -  W*  n#(X)  W±  for  all  A  6  R  . 


The  proof  of  Theorem  3.1  will  be  based  on  Theorem  2.1  and  a  series 

of  lemmas  that  relate  the  grating  propagators  A  and  A#  to  the  corresponding 

families  of  reduced  propagators  A  and  A  ,  -1/2  <  p  <  1/2. 

p  o  ,p  - 

The  Mapping  U.  As  a  first  step,  the  correspondence  Introduced  In 
[3,  (6.10)]  will  be  extended  to  a  unitary  mapping 

U  :  L2 (G)  -►  L2((-1/2,1/2],L2 (G>).  To  see  how  this  may  be  done  note  that 
if  f  €  L2  (G)  then 
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(3.5) 


f  (x+2ir4,y)  L  €  Lj(£i)  for  all  4  6  Z,  and 


(3.6) 


L  lf(*+2lT2’,‘)lL,(£i)  "  lf,L,(G)  < 


SUEZ 


Hance  Che  Plancherel  chaory  in  Che  Labesgue  space  Lj  ( (-1/2 ,1/2]  ,L2  (12) ) 
Implies  chac  Che  Fourier  series 


(3.7)  F(x,y,p)  -  l  e"2iri£p  f  (*+2rr4,y)  |Q 

46  Z 

converges  in  this  space  and  Parseval's  relaclon  is  valid.  Combining  chis 
result  and  (3.6)  gives 

(3.8)  ,F,Ll((-1/2,1/2J,L2(fi))  "  '  f '  Lj  (G) 
for  all  f6L, (G) . 

Lemma  3.2.  The  mapping  U  :  L2(G)  -  L2  ((-1/2,1/21  ,L2  (£1) )  defined 
by  Uf  ■  F  and  (3.7)  is  unicary. 

Proof.  The  preceding  discussion  Implies  chac  Uf  is  defined  for 
all  f  €  L2 (G)  and  U  is  isomecric.  The  surjeccivicy  of  U  follows  from 
Che  Plancherel  chaory.  Indeed,  every  F  €  L2((-l/2,l/2] ,L2(ft))  has  a 
Fourier  developmenc 


(3.9)  F(X,p)  -  l  e~i,rUp  F.  (X)  , 

SUEZ  * 

convergenc  in  L2 ( (-1/2, 1/2] ,L2 (fi)) .  The  Fourier  coef ficienCs in  (3.9)  are 
defined  by  che  Bochner  incegrals 
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(3.10)  FC(X)  -  f17  e27rUp  F (X,p)  dp  €  L2(G) 

^-1/2 

and  Parseval’s  relation  holds: 

(3.11)  IFI  ?  ((  .  .  -  l  IFPI.2  ,0*  . 

"2  ( C*i/2  »  y  2  ]  \U) )  ~  Lj  (**) 

Thus  to  construct  f  ■  U-1  F  one  need  only  require  that 
f(x+2iri,y)|fl  -  F£(X)  €  L2  (fl) ,  or 

(3.12)  f(X)|^(Jl)  -  FJt(x-2ir2,,y)  for  all  l  e  Z  . 

Parseval's  relation  then  guarantees  that  f  €  L2 (G)  and  (3.12)  Implies 
that  Uf  «  F. 

The  next  lemma  makes  It  possible  to  construct  operators  Y(A)  from 
the  corresponding  reduced  operators  Y (Ap)  with  p  €  (-1/2, 1/2]. 

Lemma  3.3.  For  all  bounded  Borel  functions  Yd)  defined  for 
X  >  0  and  for  all  f  €  L2(G)  one  has 

(3.13)  (U  Y(A)f) (• ,p)  -  Y(Ap)  Uf (•  ,p)  €  L2(fl) 

for  almost  every  p  e  (-1/2, 1/2]. 

Proof.  The  result  will  be  derived  from  the  R-B  wave  expansions 
for  A  and  Ap  and  the  corresponding  Plancherel  relations.  To  this  end 
let  f,g  €  L2(G)  and  write  Uf  ■  F  and  Ug  ■  G.  Moreover,  assume  that 
♦+g  *  8+  e  L2°n(R2 ) .  Then  Lemma  3.2  and  the  results  of  [5,  §5  and  §6] 
imply 


.  * 
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(U  ¥(A)f,G)  -  (¥(A)f,g)  -  (*+  ¥(A)f,  4+g)-|  2  ¥(«*(?))  ?+(P>  g+(P)dP 


(3.14) 


rnrt-1/2 

i  r 

m£Z  Jm-l/2 


¥(uz(p,q))  ?+(p,q)  g+(p,q)  dq  dp 


rl/2  _ _ 

-  I  ¥(ujz(p4in,q))  '£  (p4m,q)  g. (p+m,q)  dq  dp 

mSZ  ' -1/2  '° 


1/2 

-1/2 


r 


f°°  rTTTT— 


<¥(Ap)  F(-tp).G(*.p))  (fl)  dp  . 

Note  that  the  hypothesis  g+  €  implies  that  the  m-summation  and 

interval  of  q-integration  in  (3.14)  are  finite.  Moreover,  since  such 
functions  ^  are  dense  in  L2(R*)  the  relation  (3.14)  holds  for  all 
G€  L2  ((-1/2.1/2]  ,Lj  (0)  ) .  On  taking  G(X,p)  -  Gj(X)  G2  (p)  in  (3.14) 
where  Gt  €  L2 (0)  and  G2(p)  €  L2  (-1/2, 1/2]  are  arbitrary  one  gets  (3.13). 

The  mapping  U  obviously  depends  on  the  grating  domain  G  :  U  ■  Uq. 
In  the  special  case  G  »  R*  let  U0  ■  U^2 .  With  this  notation  one  has 
Lemma  3.4.  The  operators  U,  U0 ,  Jc  and  satisfy 

(3.15)  U„  JG  -  U  . 

Proof.  The  definition  of  U  implies  that 


(3.16)  (Jfl  Uf)(X,p)  -  e'21TUp  (f(**2ir£,y)|n)  , 
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(3.17)  (U.  Jfl  f)(X,p)  -  l  e'2ffUp  (J  f(x+2ir£,y)|  ), 

u  1€Z  o 

for  all  f  £  L2(G).  These  obviously  define  Che  same  function,  which 
implies  (3.15). 

The  next  lemma  will  be  used  to  relate  the  wave  operators  for  A 

and  A.  to  those  for  A_  and  A. 
o  P  o  ,p 

Lemma  3.5.  For  all  f  €  L2 (G)  one  has 

(3.18)  U0(*9  f)(-,p)  -  ♦*  *±  Uf (• ,p)  €  L2(B.) 

for  almost  every  p  €  (-1/2, 1/2]. 

Proof.  The  relation  [5,  (6.11)]  can  be  written 

(3.19)  (4  f)(p-Hn,q)  -  (4+  F(*,p))  (q)  «  (4,  Uf(*,p))  (q)  . 

i  ,  p  m  —  >  P  111 

The  relation  was  proved  in  [5]  for  all  f  £  L2om(G).  However  (3.19),  as 
a  relationship  in  £  ®  L2(R0),  extends  immediately  from  the  dense  set 
L2om(G)  to  all  of  L2(G).  In  particular,  specializing  (3.19)  to  G  ■  R„ 
fives 

(3-20)  (*0  f,)(p+m,q)  -  (4^p  U,  f„(* ,p))a(q) 

for  all  f0  €  L2(Rg).  Substituting  f0  ■  $+  f  in  (3.20)  gives 

(3.21)  <*0,p(uo  **  V)<*.p)>B(q)  "  V(p'Hn»q)  "  <#±, p  uf(*.p»m(q> 

in  Z  *  L2(R0),  by  (3.19).  Thus 

(3.22)  *  (U9  4*  4+f)(-,p)  -  4+  n  Uf(-.p) 

for  almost  every  p  €  (-1/2, 1/2].  (3.22)  is  equivalent  to  (3.18). 
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Lemmas  3. 2-3. 5  will  now  be  shown  Co  imply  Theorem  3.1.  The  main 
seep  in  che  proof  is  described  by 

Theorem  3.6.  For  every  h  S  L2 (G)  and  H  ■  Uh  one  has 


-itA1/2  _-itAj/2  **  *  t_R  2 


(3.23) 


r 

^  -1/2 


-itAl/2 

Jfl  6  P  H(< 


-itAl/2 

>p)-*  ‘'p  4*.P  4=.p  e(--p)Il1(b,)  dp  • 


Proof.  Lemma  3.2  implies  that 


I.  -it  A 1/2  ,  -itA^2  .  .1 2 

lJoe  h-a  *•  4±  h!L»:> 


(3.24) 


1/2 

-1/2 


MJg  e' 


-itA 


1/2 


Moreover,  Lemmas  3.4  and  3.3  imply 

-itA 1/2 

(U0  JG  e  h)(*,p) 

(3.25) 


-itA1/2 

(Jfi  U  e  1CA  h)(.,p) 
-itA1/2 

e  P  U  h(.,p)  . 


Finally,  Lemmas  3.3  and  3.5  imply 

(U„  e_itA^2  $±h)(-,p) 

(3.26) 


-itA1/2 

e  ®,P  U0($;  $±h)(.,p) 


-itA1/2 
o  .P 


$±,p  U  h(*’p)  • 


Combining  (3.24),  (3.25)  and  (3.26)  gives  (3.23). 

Proof  of  Theorem  3.1.  Lemma  3.2  implies  that  H(*,p)  €  L2 (ft)  for 
almost  every  p  €  (-1/2, 1/2].  Hence  the  integrand  on  the  right  hand  side 


( 


of  (3.23)  tends  to  zero  when  t  ■*+“  by  Theorem  2.1  (see  (2.31)). 
Moreover,  the  operators  appearing  in  the  integrand  are  all  bounded 
with  bound  1  and  hence  one  has 
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* 


H(*  ,p)  -  e 


°»P  $* 


H( 

o,p  ±,p 


•p>L<bj  i  2 


L,  (fl) 


(3.27) 


for  all  t  6  R  and  almost  every  p  €  (-1/2, 1/2].  Thus  the  existence  of  W+ 
and  W_,  and  the  relation  (3.3),  follow  from  (3.23)  and  Lebesgue's 
dominated  convergence  theorem.  The  final  statement  of  Theorem  3.1, 


equation  (3.4),  follows  from  (3.3)  and  the  eigenfunction  expansions  for 
A  and  A„ ,  exactly  as  in  the  proof  of  Theorem  2.1. 


I 


§4 .  Asymptotic  Wave  Functions  and  Energy  Distributions. 

In  this  section  the  existence  of  the  wave  operators  W+  is  shown 
to  imply  that  transient  wave  fields  in  grating  domains  G  are  asymptotically 
equal  in  the  energy  norm,  for  t  -*■  +°°,  to  transient  wave  fields  in  the 
degenerate  grating  domain  R* .  The  latter  are  then  shown  to  be  the 
restrictions  to  R2  of  free  waves  in  R2 .  Such  free  waves  possess  asymptotic 
wave  functions  in  the  sense  of  the  author's  monograph  on  scattering  by 
bounded  obstacles  [3].  These  results  are  shown  below  to  imply  that 
transient  wave  fields  u(t,X)  with  finite  energy  in  grating  domains 
possess  asymptotic  wave  functions 

(4-1)  uk(t,X)  *  r  ^2  Fk(r-t,0)  ,  k  «  0,1,2  , 


(where  X  =  (r  cos  0,r  sin  9))  such  that  if  (t,x,y)  *  (X0,Xx,X2)  and 
Dk  =  for  k  -  0,1,2  then 

(4.2)  lim  H Dj^u ( t , • )  -  u^t»*^L2(G)  *  °»  k  *  O*1*2  • 

Moreover,  the  waveforms  Fk(T,9)  are  calculated  from  the  initial  state 
f(X),  g(X)  of  u(t,X).  Finally,  (4.2)  and  the  results  of  [3,  Ch.  8]  are 
used  to  calculate  the  asymptotic  distributions  of  energy  for  transient 
wave  fields  in  grating  domains. 

The  starting  point  for  the  calculation  of  the  asymptotic  wave 
functions  (4.1)  is  the  complex-valued  wave  function  v(t,X)  defined  by 
(1.12),  (1.13).  The  existence  of  the  wave  operator  W+  defined  by  (3.1) 
implies  that 


(4.3) 


-itA1^2  -itA^2 
s  ltA  h  ~  e  ltA°  W. 


t  -►  , 
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2< 


t 

I 

| 


■j 

I 


< 


in  the  sense  of  convergence  in  L2(G).  Moreover,  if  h  €  D(A^2  )  then  the 

analogue  of  (4.3)  holds  for  the  first  derivatives.  This  result  may  be 

formulated  as  a  generalization  of  the  corresponding  result  for  exterior 

domains  [3,  Theorem  7.5],  as  follows. 

Theorem  4.1.  Let  G  satisfy  the  hypotheses  of  Theorem  3.1  and 

let  h  €  D(A1/Z  ).  Then  v(t,*)  *  e  h  is  a  solution  wFE  in  G, 

1/2 

h*  «  W+  h  €  D(Aj/2  ),  v0(t,‘)  *  e  itA°  h+  is  a  solution  wFE  in  Rg  and 

(4.4)  lim  QDk  v(t,*)  -  Dk  v0 (t, *)l^2  ^  =  0  for  k  =  0,1,2  . 

The  proof  is  precisely  the  same  as  the  one  for  exterior  domains 
given  in  [3]  and  is  therefore  omitted. 

The  initial  state  h+  *  W+  h  for  the  wave  field  v„(t,X)  satisfies 

(4.5)  *  $0  h^  =  h  *  h_ 

by  (3.3).  Thus  v0  has  the  R-B  wave  representation  (see  (1.14)) 


(4.6) 


v0(t,X) 


U>„(X,P) 


-itaj(P) 

e 


h_(P)  dP  . 


To  show  that  v0(t,X)  has  a  continuation  to  a  wave  field  wFE  in  R2  the 
Neumann  and  Dirichlet  cases  will  be  treated  separately. 

The  Neumann  Case.  Here  one  has 

(4.7)  VX,P)  "  'l'o(X.P)  -  ^  eiPx(eiciy  +  e_iqy) 

and  substitution  in  (4.6)  gives,  after  a  simple  transformation. 


(4.8) 


v0 (t,X) 


1_  [  i(xp+yq-tw(p,q)) 

27T  Jr2 


h0(p,q)  dpdq 


(; 

J 
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00 


where  Che  function  v 

(4.15)  H0(t,9)  - 


k 


are  given  by  (4.12)  with  waveforms  defined  by 

*00 

eiTU  h  (oj  cos  9,0)  sin  9)  (-iio)  dto  , 
o 


(2ir)  1/2 


convergent  in  L2(R  x  [0,ir]),  and 

(4.16)  (t, 9)  =  -H0(t,9)  cos  9  ,  H2(t,9)  =  -H0(x,9)  sin  9  . 

Equations  (4.14)  follow  from  (4.4)  and  (4.13)  by  the  triangle 
inequality.  Equations  (4.13),  (4.15)  and  (4.16)  follow  directly  from 
[3,  Ch.  2];  see  the  proof  of  [3,  Theorem  2.10]. 

To  obtain  corresponding  results  for  the  real-valued  wave  field 
u(t,X)  generated  by  the  initial  state  f,g  one  need  only  take  the  real 
part  of  v(t,X)  and  use  equation  (1.13)  which  relates  h  to  f  and  g.  This 
leads  to 

Theorem  4.3.  Let  G  satisfy  the  hypotheses  of  Theorem  3.1.  Let 
f  €  D(A^2  )  and  g  €  L2  (G)  and  define  asymptotic  wave  functions 

(4.17)  u*(t,X)  -  r"1/2  Fk(r-t,9)  ,  k  -  0,1,2  , 
by 


F0(T,8) 


(4.18) 


Re 


(2ir) 


T TT 


iTU) 


1  0 


[g_(w  cos  9,w  sin  9)  -  iwf_((D  cos  8,w  sin  9)  ]  (-iw)  du)j  , 


convergent  in  L2(R  x  [0,ir]),  and 


(4.19) 


Fi(t,0)  -  -F0(t,9)  cos  8  ,  F2(t,8)  -  -Fj (t , 8)  sin  8)  . 
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Then  the  solution  wFE  (1.10)  generated  by  f  and  g  satisfies 


(4.20)  lim  lDk  u(t,*)  -  uk(t,-)lL  (G)  -  0,  k  -  0,1,2  . 

Proof.  To  begin,  assume  that  g  s  D(A  172  )  and  define  h  by 

(1.13).  Then  h  e  D(Al72  )  and  Corollary  4.2  is  applicable  to 
-itA172 

v(t,*)  -  e  h.  Moreover, 


(4.21) 


|P|  h_  (P)  -  |P|  f_(P)  +  i  g_  (P) 


which  implies  that  *  Re  Hk  and  hence 


(4.22) 


uk(t,X)  -  Re  (vk(t,X)}  ,  k  -  0,1,2  . 


Thus  (4.20)  follows  from  (4.14)  by  the  triangle  inequality.  To  remove 
the  restriction  that  g  £  D(A  l7z  )  note  that  D(A  1/,z  )  is  dense  in  L2(G), 
by  the  spectral  theorem.  Moreover,  one  has  (see  [3,  Theorem  2.5]) 


i  "''"ct. -)l1_z (R,)  i  IH„I  -  IH0I 


(4.23) 


|  ^  |o  cu 3  |h_(uj  cos  0,0)  sin  6)  | 2  d0  do>| 
ff  -  U/2 

y  ipi2  ih.(p>i2  dpj  -i|pi  h.(P)iLi(R2) 


<  I  |p  | f_ (P)S  +  IgJ  -  1A1/2  fl  +  Igl  . 


Similarly, 


(4.24) 


Iuk(t,*)lL  <  lu"(t,*)l  <  SAl/2  fl  +  Igl  ,  k  -  1,2 


Finally,  the  conservation  of  energy  theorem  implies  that 
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(4.25) 


5Dku(t,‘)llL2(G)  <  E(u,G,t)1/2  -E(u,G,0)l/2 

-  (I  A 1/2  fi2  +  9  gl  2> 1/2  <  lAl/2  fl  +  Igl  ,  k  -  0,1,2  . 


It  follows  from  (4.23),  (4.24)  and  (4.25)  that  (4.20)  can  be  extended  to 
all  f  6  D(A^2  )  and  g  6  L2(G)  by  a  well-known  density  argument  (see,  e.g., 
(3,  Proof  of  Theorem  2.6]). 

Theorem  4.3  permits  the  extension  to  grating  domains  of  the 
results  on  asymptotic  energy  distributions  in  exterior  domains  given  in 
[3].  The  principal  results  are  formulated  below.  The  proofs  are 
identical  to  those  of  [3]  and  are  therefore  omitted. 

Corollary  4.4  (Scattering  into  Cones).  Let 

(4.26)  T  »  (X  -  (r  cos  0,r  sin  0)  :  r  >  0  and  0  e  r0) 


where  r0  is  a  Lebesgue-measurable  subset  of  [0,tt],  and  let  X0  e  R2.  Then 
under  the  hypotheses  of  Theorem  4.3  the  limit 


(4.27)  E°°(u,G  n  (T  +  X„))  -  lim  E(u,G  n  (T  +  X0),t) 

t-*+» 


exists  and 


(4.28)  E*(u,G  n  (r  +  X0))  =  ||P|  f  (P)  +  i  g  (P)  l 2  dP  . 

jr 

Corollary  4.5  (Transiency  of  Energy  in  Slabs).  Let 

(4.29)  Z  -  {X  :  dt  <  X  •  X0  <  d2} 

where  d}  and  d2  (>  dj)  are  constants  and  X„  e  R2  is  a  unit  vector.  Then 
under  the  hypotheses  of  Theorem  4.3  one  has 
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(4.30)  lim  E(u,G  O  E,t)  -  0  . 

Note  that  Corollary  4.5  implies  the  transiency  of  the  energy  in 
bounded  sets  since  every  bounded  set  K  C  R2  is  contained  in  a  slab 
(4.29). 


/ 


15 .  Construction  and  Structure  of  the  S-Matrix. 

The  scattering  operator  associated  with  the  pair  A,  A0  is  the 
linear  operator  S  :  L2(Rj)  ■*  L2(Rg)  defined  by 

(5.1)  S  -  W+  W*  . 

The  corresponding  operator  in  L2(R*)  defined  by 

(5.2)  S  -  S  ** 


is  the  Heisenberg  operator,  or  S-matrix,  for  the  pair  A,  A0 .  From  the 
representation  W+  -  $_  of  Theorem  3.1  one  has 

(5.3)  S  -  4>_  . 

The  unitar ity  of  W+  and  <60  imply  that  S  and  S  are  unitary  operators  in 
L2(Rq).  The  purpose  of  this  section  is  to  calculate  S.  Specifically, 

A 

it  will  be  shown  how  S  can  be  constructed  from  the  scattering  coeffi¬ 
cients  (c^(p,q)}  of  the  R-B  waves  ^+(x,y,p,q)  and  the  relationships 
among  these  coefficients  imposed  by  the  unitarity  of  S  will  be  determined. 
The  role  of  the  S-matrix  in  the  scattering  of  transient  fields  by 
gratings  will  be  developed  in  §6. 

If  h  €  L2(G)  then  (5.3)  implies  that  the  functions  h+  ■  *+  h  and 

/v 

h_  »  <t>_  h  satisfy 


(5.4) 


h_  -  S  h+  . 


A  ^ 

Thus  S  may  be  calculated  by  calculating  the  relationship  between  h_  and 

/v 

hj_.  This  will  be  done  by  using  the  incoming  and  outgoing  R-B  wave 


31 


32 


-itAI/2 

representations  of  v(t,')  ■  e  h  to  calculate  in  two  different  ways 

the  asymptotic  wave  function  in  L2(G)  associated  with  v(t,*);  say 

(5.5)  v°°(t,X)  -  r  l^2  H(r-t, 0)  . 


The  function  H  S  L2(R  x  [0,tt])  is  uniquely  determined  by  the  condition 


(5.6) 


lim  I v( t , • )  -  v  (c , • ) »  -  0  ; 

^-►+00  “2  'v'' 


see  [3,  Theorem  2.5].  The  equality  of  the  representations  of  H  obtained 
from  the  incoming  and  outgoing  representations  of  v(t,*)  provides  the 

A  A 

required  relationship  between  h_  and  h+. 

First  Calculation  of  H.  Theorem  3.1  implies  that 


(5.7)  lim  lv(t,*)  -  v0(t,*)lL  (r)  ■  0 

oo 

where  v0(t,*)  ■  e  hg  Is  the  wave  function  in  L2(R*)  of  Theorem  4.1. 

Proceeding  as  in  the  proof  of  Corollary  4.2  one  shows  that  (5.5),  (5.6) 
hold  with 


(5.8)  H(t,  0)  -  ^2tt)  i77~  ]„  eiTai  h_(w  cos  0,u>  sin  0)  (-iui) l^2  dw  . 

The  convergence  v0(t,*)  -  v“(t,»)  ■+■  0  in  L2(R2)  was  proved  in 
[3,  Theorem  2.6]. 

A  Classification  of  the  R-B  Waves.  The  second  calculation  of  H 
will  be  based  on  the  outgoing  representation 

(5.9)  v(t,X)  -  l.i.m.  f  if;,(X,P)  e"itw(P)  h,(P)  dP  . 

Jr2  +  + 

The  R-B  wave  <J>+  has  the  expansion  for  y  >  h  [5,  (2.26)] 
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iCp^x+q^y) 


(5.10) 


^(x.y.p.q)  -  (2tt)~1  ei(px”qy)  +  £  ct(p,q)  e 

(fH-Jl)2<p2+q2 


(p+l)2>p2+q2 


where 


(5.11)  (p^.q^)  *  (p  +  £.{p2  +  q2  ~  (p  +  £)2}1/f2 )  . 

The  first  sum  in  (5.10)  is  a  superposition  of  a  finite  number  of  outgoing 
plane  waves,  while  the  second  sum  is  an  exponentially  decreasing  function 
of  y  for  (x,y)  €  Rg  -  E  [6,  Lemma  9.3].  In  the  calculation  of  the 
asymptotic  wave  function  (5.5)  from  (5.9)  and  (5.10)  a  difficulty  arises 
because  the  number  of  terms  in  the  first  sum  varies  with  (p,q)  6  . 

This  number  changes  at  the  points  (p,q)  6  E  and  is  constant  on  the 
components  of  the  set  R2  -  E.  It  will  therefore  be  convenient  to  classify 
the  R-B  waves  by  means  of  these  components.  Note  that  (p,q)  €  E  if  and 
only  if  q  >  0  and 

(5.12)  q2  =  p2  +  q2  -  (p  +  l)2  -  0  ,  l  €  Z  -  {0}  . 

The  set  so  defined  is  the  portion  lying  in  R2  of  the  parabola  with 
focus  at  (0,0)  and  vertex  at  (-Jt/2,0).  The  curves  and  it  are  disjoint 
if  >  0  and  intersect  orthogonally  if  Jim  <  0.  Thus  if 

"  Ro  n  t(p’q)  :  |P  +  “I  <  /p2  +  q2  <  |p  +  m  +  1 | >  , 

0_n  "  Ro  n  ^  <P»<J>  ’■  |p  -  n|  <  /p2  +  q2  <  |p  - 


(5.13) 


n  -  1 1 }  , 
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where  m,n  ■  0,1,2, ••*,  then  0  is  the  domain  between  tt  and  tt  a  ,  0  is 

m  m  nrr  1  — n 

the  domain  between  tt  and  tt  and  the  sets 

-n  -n- i 


(5.14) 

are  the  components  of  R*  -  E  : 


0  -  0  n  0 

m,n  m  -n 


m,n  -  0,1,2, •••, 


(5.15) 


K-e 


0 


—  n  m,  n 
m,n-0 

Note  that  (p,q)  60  if  and  only  if  the  expansion  (5.10)  of  <p.(x,y,p,q) 
m  y  n  t 

contains  exactly  m  +  n  +  1  outgoing  plane  waves  with  the  propagation 
directions  (P£»<l£)>  -n  <  l  <  m. 

Second  Calculation  of  H.  In  calculating  S  it  will  suffice  to 


A  A  A  _  A 

determine  S  h+  for  functions  h+  of  a  dense  set  in  L2(R*)  because  S  is 
known  to  be  unitary.  For  this  purpose  it  will  be  convenient  to  use 
functions  h+  €  C0(RQ  -  c) .  For  such  functions,  supp  h+  is  a  compact 
subset  of  the  set  (5.15).  Hence,  supp  h+  meets  only  finitely  many  of 
m  n  and  each  component  of  supp  h+  lies  in  one  of  these  sets. 

A  A 

Thus  in  calculating  S  h+  it  will  be  enough  to  consider  the  case  where 


(5.16)  supp  h ,  ■  K  C  0  ,  m  and  n  fixed  . 

t  in,  n 

The  case  of  a  general  h+  €  c"(R*  -  E)  may  then  be  obtained  by  super¬ 
position.  With  this  hypothesis  the  wave  function  (5.9)  becomes 


(5.17) 


v(t,X) 


KCX.P)  e-ic“<P> 
K 


h+(P)  dP  . 


The  asymptotic  wave  function  (5.5)  for  v(t,X)  will  be  calculated 
from  (5.17)  by  substituting  the  expansion  (5.10)  and  determining  the 


35 


behavior  for  c  ■»  +*  of  the  terms  in  the  resulting  sum.  For  this 
purpose  a  bound  is  needed  for  the  remainder  in  (5.10)  that  is  uniform 
in  (p,q)  €  K.  Thus  a  refinement  of  [6,  Lemma  9.3]  is  needed  since  the 
latter  is  valid  for  fixed  p  only.  The  following  generalization  of 
[6,  Lemma  9.3]  will  be  proved. 

Lemma  5.1.  Define  the  remainder  a+(X,p,q)  for  all  Xe  G  and 

(p.q)  €  0  _  by 

m,  a 


(5.18)  Mx.P.q) 


(  2tt) 


gi(px+qy) +  r 


cj (p,q) 


i(P»x±Qoy) 

e  +  a±(X,p,q) 


Then  for  each  compact  set  K  C  0  and  each  r'  >  r  >  h  there  exist 

m,n 

constants  U  -  U(K)  >0  and  C  -  C(K,r',r)  such  that 

(5.19)  Ja±(X,p,q) |  <  C  e“Uy  for  all  X  6  R* ,  and  (p,q)  6  K  . 

Proof .  Only  the  case  of  a+  will  be  discussed  since  the  other 
case  then  follows  from  the  relation  [5,  (2.25)].  The  proof  will  parallel 
that  of  [6,  Lemma  9.3].  Note  that  [5,  (6.5)]  implies  that 

(5.20)  0+(x,y,p,q)  -  e27TiZP  p+(x  _  2ir£,y,p,q)  ,  (x,y)  6  , 

where  p+  is  defined  by  [6,  (9.39)]  with 

(5.21)  L'  ■  {£  €  Z  :  i  <  -n  -  1  or  H  >  u  +  1} 

for  all  (p,q)  e  0^  q.  Thus  to  prove  (5.19)  it  is  enough  to  show  that 

(5.22)  |p+(X,p,q)j  <  C  e_My  for  all  XS  fl  ,  and  (p,q)  €  K  . 
Proceeding  as  in  [6,  Proof  of  Lemma  9.3],  one  has 


36 


(5.23)  k;A(r,p,q)  |2  <  C2(2iT)"l]<{)|(*,p,q)I  ^<h  ^  for  all  i  €  L' 

where  C0  *  C0(h,r).  Now  the  right-hand  side  of  (5.23)  is  a  continuous 
function  of  (p,q)  SR*-  E  by  [5,  Theorem  6.1].  Thus  there  exists  a 
Cj  -  Cj(K,r)  such  chat 

(5.24)  |<t>|£(r,p,q)  |  <  Cl  for  all  (p,q)  €  K  and  l  €  L'  . 

Next,  since  K  is  a  compact  subset  of  0  there  exist  constants 

m,  n 

U+  -  M+(K)  >  0  and  M_  -  U_(K)  >  0  such  that 

(p  +  Jl)  2  -  p2  -  q2  >  U*  for  all  (p,q)  €  K  and  l  >  m  +  1  , 

(5.25) 

(p  +  £) 2  -  p2  -  q2  >  u2  for  all  (p,q)  6  K  and  l  <  -n  -  1  . 


whence 


(5.26)  {(p  +  £) 2  -  p2  -  q2}1/2  >  y(K)  -  Min  (p+(K) ,y_(K))  >  0 

for  all  (p, q)  €  K  and  l  €  L'.  It  follows  that  for  all  X  e  Q  ,  and 

r 

(p,q)  S  K  one  has 


(p  (X,p,q) |  <  l  (y,p,q)| 

££L' 

(5.27) 

-  I  U+o<r»P»q)  I  exp  (-(y-r)((p+^)2-p2-q2)1/2} 

jt<=L ' 

<  Ct  l  exp  {-(y-r) ( (p+£) 2-p2-q2) l/2} 

2€L' 
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t 


<  C,  l  exp  (-(y-r ') ((p+2) 2-p2-q2) 1/2}  exp  {-(r '-r)  ((p+fc)  2-p2-q2) l/2} 
(5.27  cone.) 

<  C.  e-(y-r’)U(K)  I  exp  {-(r’-r)((p+Jl)2-p2-q2)l/2}  . 

2£L' 


Now 


(5.28)  Z(r'-r,p,q)  *  £  exp  {-(r ’-r)  ((p+Z)  2-p2-q2)  l/2} 

JteL’ 

is  a  continuous  function  of  (p,q)  £  0  and  hence  for  each  compact 

0,0 

K  C  0  there  is  a  constant  M(r’  -  r,K)  such  that 
m,  n 

(5.29)  Z(r'-r,p,q)  <  M(r'-r,K)  for  all  (p.q)  £  K  . 

Combining  (5.27),  (5.28)  and  (5.29)  gives  (5.22)  with 

(5.30)  C  -  Cj(K,r)  er’U(K)  M(r'-r,K)  . 

Second  Calculation  of  H  (continued).  Substitution  of  (5.18) 
into  (5.17)  gives  the  decomposition 

.  m 

(5.31)  v(t,X)  -  vln(t,X)+  l  v°uC(t,X)  +  va(t,X)  ,  t  6  R  ,  X€  G  , 

£«-n 

where 


(5.32) 


(5.33) 


out 


vin(t,X) 
(t,X)  -  f 


K 


-  ±  j  ,1<P*-«I7-Mp.q»  h+(p,q)  dpdq  , 
R 

i(p,x+q,y-ta)(p,q)) 

e  ca(p,q)  h+(p,q)  dpdq 


and 
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(5.34)  v0(t,X)  “  |  e“it(JJ(P» q)  0+(X,p,q)  h+(p,q)  dpdq  . 

K 

Recall  Chat  by  assumption  h+  £  C0(Rg  ~  E)  satisfies  (5.16)  and  c^(p,q) 

£  C(Rq  -  E)  (see  (2.23),  (2.24)).  Thus  the  integrands  in  the  above 
integrals  are  all  continuous.  The  second  calculation  of  H  will  now  be 
carried  out  by  calculating  the  asymptotic  wave  function  in  L2(G)  of  each 
term  on  the  right-hand  side  of  (5.31). 

The  Partial  Wave  vin(t,X) .  The  change  of  variables  (p',q') 

■  (p,-q)  in  (5.32)  gives 


(5.35) 


vin(t,X) 


_1_  (  i(p'x+q,y-tu)(p',q')) 

2ir  JK- 


h+(p\-q’)  dp'dq' 


where  K'  »  ((p',q')  :  (p,q)  €  K}  C  R2  -  R2.  Thus  vin(t,X)  is  a  free 
wave  in  R2  and  hence  has  an  asymptotic  wave  function  r_1/^2  H^n(r  -  t,0) 
with  waveform  defined  by  (3,  Theorem  2.6] 


(5.36)  Hin(T,0) 


(  2-tt) 


iTW 


e  h+(u)  cos  0 ,  —oo  sin  0)  (-iu>) 


1/2 


du) 


In  particular,  Hin(T,0)  =  0  for  0  <  0  <  tt  because  Kr  «  supp  h.(p,-q) 

“  “  *T* 

C  R2  -  Rj  . 

The  Partial  Waves  v£ut(t,X).  To  interpret  these  terms  let  l  €  Z 
and  consider  the  mapping  X^  defined  by 

(5.37)  (p^.q^)  *  xg,(P*^  "  (P  +  Mp2  +  q2  -  (P  +  )  . 


X^  is  analytic  on  the  domain 
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(5.38)  D(X^)  =  { (p, q)  :  /F“+“5r  >  |p  +  £|  ,  q  >  0} 


and  maps  it  bijectively  onto  the  range 


(5.39)  R(Xa)  -  Up^.qj)  :  ^  +  >  i ~  »  9^  >  0} 


Moreover, 

(5.40) 


.-i 


t  -  x_z  ,  i  e  z  , 


and  the  Jacobian  of  X^  is 


(5.41) 


3(P£»q£)  q 
3(p»q)  q£ 


Note  that  oj(p,q)  is  invariant  under  X^  : 


(5.42) 


in 


(p^.q^)  “  ?l  +  =  p2  +  q2  *  «(p»0 


It  is  easy  to  verify  that 


(5.43) 


X£  °m,n  “  -n  <  £  <  m 


Hence  the  hypothesis  K  C  0^  ^  implies  that 


(5.44) 


x£  K  1  H  C  V*,«+!  for  -n  <  1  <  m 


In  particular,  one  has 


(5.45) 


n  »  4>  for  j  f  i 


On  making  the  change  of  variables  (p,q)  -►  (Pg^q^)  *  X^ 


(p,q)  in 


(5.33)  one  finds  the  representation 
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(5.46) 


out... 

V0  (t,X) 


X 


Kxp^yq^aXp^.q^Jj  +  q^ 

i  |ci(p,q)h+(p,q)  — | 


dpz  dq^ 


where  (p,q)  =  X_ ^ (p^, q^)  in  the  integrand.  Thus  v£ut(t,X)  is  also  a 
free  wave  in  R2  and  has  an  asymptotic  wave  function  r  l^z  H^(r-t,0)  with 
waveform  defined  by 


i*00 


(5.47) 


H0(T,6) 


(2TT)  V2 


£ixw  h^(oj,0)  dw 


and 


(5.48)  Ha(w,0)  -  27r(-iu))1/z  —  c*(p,q)  h+(p,q) 

q 


(p> q)=X_£ (to  cos  0,0)  sin  6) 


A  simple  calculation  shows  that 


(5.49) 


I  H 


i  L2(R*[0,1T]) 


(2tt)  ‘ 


ct(p,q)  h. 


q0 

.(P,q) I2  —  dp  dq 


where  q^  is  defined  by  (5.37).  In  particular,  H£  e  L2 (R  *  [ 0 , tt ] ) 
because  the  integrand  in  (5  .49)  is  continuous  on  K. 

The  Partial  Wave  va(t,X).  Equation  (5.31)  may  be  written 

(5.50)  v  (t,X)  *  v(c>x>  ~  vln(t,X)  -  l  v°ut(t,X) 

£--n 


for  all  t  €  R  and  X  €  G.  Moreover,  it  has  been  shown  that 

v(t,X)  -  r"1/z  H(r-t, 0)  +  o(l)  , 
vin(t,X)  -  o(l)  , 
v°ut(c,X)  -  r_1/z  H^r-t,©)  +  o(l) 


(5.51) 
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where  each  term  o(l)  €  L2(G)  for  all  t  S  R  and  tends  to  zero  in  L2(G) 
when  t  -*•  +°°.  These  results  imply  that 

(5.52)  va(t,X)  «  r_l/z  Ha(r-t,0)  +  o(l) 
where  o(l)  +  0  in  L2(G)  when  t  -*•  +00  and 

m 

(5.53)  Ha(T,8)  =  H(t,  9)  -  l  H£(t,8)  in  L2  (R  *  [O.ir])  . 

£=•- n 


On  the  other  hand  (5.34)  and  Lemma  5.1  imply  that 


(5.54) 


|va(t,X)|  <  Ca  e"yy  for  all  t  €  R  and  X  €  Rz, 


where 


(5.55) 


Ca  *  c 


K 


|h+(p,q) |  dp  dq 


and  y  =  y(K)  >  0  and  C  =  C(K,r',r)  are  the  constants  of  the  lemma.  The 
second  calculation  of  H  will  now  be  completed  by  showing  that  (5.52)  and 

(5.54)  imply 

Theorem  5.2.  Ho(t,9)  =  0  and  hence 

m 

(5.56)  H(t , 9)  =  l  H  (t,0)  in  L2(R  x  [ 0 , tt ] )  . 

I— n  X 

Proof .  Let  £  be  an  arbitrary  number  in  the  interval  0  <  e  <  it/ 2 
and  consider  the  sector 

T£  *  {(x,y)  -  (r  cos  0,r  sin  8)  :  £  <  0  <  it  -  e}  . 


(5.57) 


By  (5.52),  the  local  decay  of  asymptotic  wave  functions  [3,  p.  32]  and 
the  triangle  inequality,  one  has 


'Gnr, 


|va(t,X)|z  dX 


|Ha(r-t,0) |2  r-1  dX  +  o(l) 


(5.58) 


e 

®  p-e 
o 


(•tt-e 

|Ha(r-t,6) |2  d6  dr  +  o(l) 

'  £ 
r7T-£ 

=  |Ha(T,8)|2  de  dT  +  o(l) 

J-t  £  ° 


where  o(l)  **■  0  when  t  -*■  00 .  Thus  passage  to  the  limit  in  (5.58)  gives 


(5.59)  lim 

fN-°° 


Gnr, 


|va(t,X) 1 2  dX 


ir-e; 


|Ha(T,e) |2  de  dt  . 


On  the  other  hand,  writing  R2  *  {(x,y)  :  x  €  R,  a  <  y  <  b}  ,  one  has 

dy  D 


iva(t,X)|2  dX 


Gnr, 


=|  |va(t,X) |2  dX  +  |vo(t,X)|2dX 


(5.60) 


Gnr  nR2  , 
e  o,k 


T  nR2 


<  c* 


-  a 


-jyr  sin  0  ,  ,  ... 

e  r  dr  d0  +  o(l) 


r  np2 


K 


for  every  fixed  k  >  r'  >  r  by  (5.54)  and  the  local  decay  property  for  v 
Passage  to  the  limit  t  -H»  in  (5.60)  gives  by  (5.59), 


(5.61) 


fir-e 

|  |ho(t,0) |  d0  dT  <  C2 


«°  rTT-e 
-00  j  E 


^-2yr  sin  0  ,  ,, 

e  r  dr  d0 


r  np2 


% 


for  every  k  >  r'.  Note  that  C  ,  y  and  the  left-hand  side  of  (5.61)  are 


independent  of  k.  Now,  sin  0  >  sin  e  >  0  for  e  <  0  <  it  -  e  and  hence 
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-2yr  sin  8  ,  ..  . 

e  K  r  dr  d0  < 


(5.62) 


r  hr* 

£  k 


e-2yrsine  r  dr  d0 


r 


■H 


f7r’£  e-^sine  r  d6  dr  ,  (ir  .  2e) 

£  ; 


-2yr  sin  £  , 

e  r  dr 


But  the  last  integral  tends  to  zero  when  k  -*■  <*>  with  e  fixed.  Thus  (5.61) 
implies  that  Ha(t, 9)  =  0  in  R  x  [£,u-£]  and  (5.56)  follows  since 
£  €  [0,tt/2]  is  arbitrary. 

Corollary  5.3.  For  all  h  6  L2(G)  such  that  supp  h+  C  n 

=  closure  of  0  one  has  the  two  relations 
m,n 

(5.63)  h  (p,q)  =  2H  £  Cj(X_{(p,q))  h  (X  f(p,q)) 

Sb—n  *  *  ~  1  q-l 


for  almost  every  (p,q)  €  Rg  where  (p_£,q_g)  *  X_^(p,q). 

A  - 

Proof.  The  case  where  supp  h,  C  0  is  considered  first.  In 
-  +  m,n 

this  case  it  will  suffice  to  prove  (5.63)  for  functions  h  €  L2(G)  such 

*  oo  - 

that  h,  e  C.(R.  -  £)  and  supp  h ,  *  K  C  ()  since  such  functions  are 
-r  u  «  -r  m,n 

dense  in  the  subspace  of  L2(G)  defined  by  supp  h+  C  n>  For  such 
functions  the  relation  (5.56)  and  the  Fourier  representations  (5.8)  for 
H(X,8)  and  (5.47),  (5.48)  for  H^(t,9)  imply  that 

h_(o)  cos  9,0)  sin  9) 

(5.64) 

m 

cJ^(X_Ji(o)cos  0,0)  sin  0))h+(X_^(o)cos  8,o)  sin  8))  |^2_(aJc0s  Q_ft)'2}  1/2 
for  almost  every  (o),9)  £  R,  x  Making  the  substitutions 
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p  ■  a)  cos  9,  q  **  w  sin  9  in  (5.64)  gives  (5.63)  in  the  case 


supp  h+  C 


m,n 


The  second  case  of  (5.63)  can  be  derived  by  calculating  the 

-itAl/2 

asymptotic  wave  functions  for  v(t,*)  ■  e  h  when  t  •*  using 

the  method  given  above.  A  simpler  derivation  may  be  based  on  the 
relations  4<_(X,p,q)  -  ip+(X,-p,q)  and  c^(p,q)  -  c+^-p.q)  of  [5,  (2.25) 
and  (2.29)].  Indeed,  if  supp  c  n  and  S(X)  *  h(X)  then  these 
relations  imply  that  g+(p,q)  *  (p,q)  and  hence  relation  (5.63)  with 

the  upper  sign  for  g  implies  (5.63)  with  the  lower  sign  for  h. 

A 

The  Structure  of  S.  It  will  be  convenient  to  use  the  notation 


(5.65)  «.,„»>  *  Vn<P>  8<P) 

where  x  is  the  characteristic  function  of  the  set  0  .  Clearly,  the 

m,  n  m,  n 

operator  Pm  q  in  Lz(Rg)  defined  by 


(5.66) 


P 

m,n 


g 


8_  _  9 

in,  n 


m,a 


o,  1,  2,' 


is  an  orthogonal  projection  and  different  operators  of  the  family  have 
orthogonal  ranges.  Moreover,  the  relation  (5.15)  implies  that  the 
family  is  complete  because  e  is  a  null  set;  i.e.. 


(5.67) 


00 


l 

m,n*0 


P 

m,  n 


1  . 


It  follows  that  for  all  g  e  LZ(R2)  one  has 


(5.68) 


S  g 


m,n*0 


S (g_  _) 
in,  n 


Thus  S  is  completely  determined  by 
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Theorem  5.4.  For  all  g  €  L2(R?)  one  has 


m 


(5.69) 


S^m,n‘>  .  ^  ^  ®m,n^m-£,rH-£  *  an<i 

£--n 


<S-70>  <W»-t.W-)l<P>‘,)  '  2”  tl  4<X-t(l>'',)> 


Similarly,  one  has 


(5.71) 


"  J  (§*  8m,n>m-£,rH-£  ’  30,1 
£*-n 


(5.72)  (5*  sa,>n)m.t,0+i(p.q)  •  2*  ^  Htt.t<p.q)>  em,n«-l(p.t» 


In  particular,  if  supp  g  c  n  then 

a  a  m 

(5.73)  supp  S  g  U  supp  S*  C  |J  . 

£«-n 

Proof .  Equations  (5.69),  (5.70),  (5.71)  and  (5.72)  follow 
immediately  from  Corollary  5.3,  the  relations  *  S  h+,  h+  *  S*  h_  and 

A  —  ^ 

the  observation  that  when  supp  h,  C  0  then  the  £  term  in  the  sum  in 
(5.63)  has  its  support  in  (5.73)  follows  from  (5.69)  and 

(5.71). 

The  unitarity  of  S  and  (5.70),  (5.72)  impose  restrictions  on  the 
scattering  coefficients  c^.  To  calculate  them  it  will  be  convenient  to 
calculate  S*  and  S  directly  from  (5.70)  and  (5.72),  respectively.  This 

A  •$(  A 

gives  the  following  alternative  representations  of  S  and  S. 


♦ 
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Theorem  5.5.  For  all  g  €  L2(Rj)  one  has 


(5-74)  (S  *«,n)«-£,iH-i(p*q)  "  211  c-i(p’q)  gm,n(X-£(p’q)) 


and  similarly 


(5‘75)  (§  gm,n>m-£,irt-£(P‘q)  "  2Tr  cIi<p’q)  gm,n(X-£(p’ q) > 


(5.76) 


Proof.  For  all  f,g  £  L2(R*)  one  has 


(f,<s*  g)m  n)  -  (f  n,s*  g)  -  (S(f  >,g) 

m,n  m,n  m,n 


m 


"  A  (S(fm,n),gm-Jl,iH-£> 
J6“-n 


27r  l  {  cj^X-2?  fm,n^X-£^  gm-£, 
£«-n  ’  ’ 


n+£  q 


_S_ 


dp  dq 


-l 


m-£,n+£ 

by  (5.70).  On  making  the  change  of  variables  (p',q')  *  (p_^»q_Jl) 

*  X_^(p,q)  in  the  last  integral  and  noting  that  q/ q_^  ■  3(p' ,q')/9(p»q) , 


one  has 


m 


<f,(S  «>„,„>  -2-  I 

X/*~n 


c£(p,'q,)fm,n(P’*q,)  gm-£,rH-Jl(X£(p,’q,))dp,dq, 


(5.77) 


m,n 


Ro 


f(p,q) 


f  n  — - 

2lT  JLn  ^(P,q)  V£,n+£  (X^Cp.q)) ^ 


dp  dq 


because  supp  gn_^  ^^(X^)  c  n-  Since  f  €  L2(R*)  is  arbitrary,  (5.77) 
implies  that 


(S 


g)m,n<p*q)  "  2lT  l  ct(p’q)  gm-£,n+£(VP’q))  * 

n 


(5.78) 
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f 


To  derive  (5.74)  note  that  for  all  l,  m,  n,  m,  n  >  0 


(5.79) 


n^m-£,n+£,  *  ^m-£,m  ^n+£,n  ®m,  n 


where  <$.,  is  the  Kronecker  symbol.  Noting  that  6  „  —  6  . „  — 

jk  jo  n+Jt,n 

“  5 2.  m  m  5  2  n  n  "  m  n  n  5  2.  m  (5.78)  and  (5.79)  imply 
x,m-m  x.,n-n  m-m,n-n  *>,m-m 


(S*  g-  -)  (p. q)  -  2ir  6  -  -  l  c+(p,q)  6.  -  g-  -(X,  (p,q) ) 

iq  9  q  in  5  n  iq  ni  ^  n—ii  «  x»  m~oi  ni«ri 

X-- - n 


(5.80) 


2tt  6 - 


m-m,n-n  Cm-m^P,C^  ®m,n^m-m^p,<^ 


This  clearly  is  zero  unless  m-m»n-n»2,  where  -n  <  l  <  m,  which 
implies  (5.71).  Moreover,  setting  (m,n)  *  (m  -  £,n  +  £,)  in  (5.80)  gives 
(5.74).  The  proof  of  (5.75)  is  obtained  by  the  same  method,  beginning 
with  (5.72). 

The  two  representations  of  S  and  S*  of  Theorems  5.4  and  5.5  hold 
for  arbitrary  g  €  L2 (R*) .  It  follows  that  the  scattering  coefficients 
must  satisfy  the  relations 


(5-81)  q  c^(X_£(p,q))  »  q_£  cf£(p,q) 

A 

for  all  (p, q)  €  ^ Moreover,  the  unitarity  of  S  and  Theorems  5.4 

and  5.5  imply 

Theorem  5.6.  The  scattering  coefficients  c^  of  the  R-B  waves 
'l'±(X,p,q)  satisfy  the  identities 


(5.82) 


2  +  +  -2 
l  C£(p,q)  c-_k  (Xk(p,q))q£  -  (2ir)  q  6k  #  ,  and 
£*-n  ’ 


i 
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(5.83) 


+  +  —  1 

l  (QCX^p.q))  c“_4(XJl(p,q))q^ 

£»-n 


(27r)_2q~1 


for  all  (p,q)  £  0  and  all  k  such  that  -n  <  k  <  m. 
r  m,n 

These  properties  may  be  verified  by  simple  calculations  using 
the  relations 


(5.84) 


(S(fm,n)’S(8m-k,trt-k))  *  6k, « ^m.n*^,^  *  and 


(5.85) 


o^m,n’®m,n^ 


and  the  constructions  of  S  and  S  described  in  Theorems  5.4  and  5.5. 
Relation  (5.83)  also  follows  from  relations  (5.81)  and  (5.82). 

It  is  well  known  in  the  theories  of  scattering  by  potentials  and 

A 

by  bounded  obstacles  that  the  S-matrix  S  is  a  direct  integral  of  a 

A 

family  of  unitary  operators  S(u>)  that  act  on  the  "energy  shell" 
p2  +  q2  ■  oj2.  The  analogous  property  of  the  S-matrices  for  diffraction 
gratings  is  evident  from  Theorem  5.4  and  the  properties  of  the  mappings 
X^.  The  operator  S(w)  in  this  case  is  given  by  (cf.  (5.75)) 


A  1  _ 

(5.86)  S(oj)  g(u)  cos  9,u  sin  9)  *  2tt  Y  c0  (w cos  8,u)  sin  9)  g(X„  (u>  cos  8,u>  sin  9) ) 

fc“-m  1 


when  supp  g  C  () 

r  °  m,n 

tion  with  supp  s  C  (9 
written 


If  s(9)  *  g(u>  cos  9, w  sin  9) 

:  (uj  cos  8,u>  sin  8)  g  0  } 

m,  n 


is  an  arbitrary  func- 
then  (5.86)  can  be 


A  H 

(S(uj)s)(9)  -  2tt  £  c^Ou  cos  8, to  sin  6)  s<0^> 
2.*-m 


(5.87) 
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where  8^  ■  9^(uj,0)  is  defined  as  the  unique  angle  such  that  0  <  9^  <  it  and 

(5.88)  X^(o>  cos  9,oi  sin  9)  •  (u  cos  0^,0)  sin  9^)  . 

For  general  s  £  L2(0,it),  S(u)s  is  obtained  from  (5.87),  (5.88)  by  super¬ 
position.  The  unitarity  of  S(w)  in  L2(0,tt)  can  be  verified  by  direct 
calculation  using  (5.87),  the  analogue  for  S*(o»)  and  Theorem  5.6. 


§6 .  The  Scattering  of  Signals  by  Diffraction  Gratings. 

The  results  of  §4  and  §5  are  applicable  to  the  echoes  that  are 
produced  when  signals  generated  by  localized  sources  are  scattered  by  a 
diffraction  grating.  The  structure  of  such  echoes  is  analyzed  in  this 
section.  Most  of  the  section  deals  with  the  case,  often  realized  in 
applications,  where  the  sources  are  far  from  the  grating.  In  particular, 
it  is  shown  that  in  this  case  the  influence  of  the  grating  on  the  echoes 
is  completely  described  by  the  S-matrix. 

It  will  be  assumed  that  the  sources  of  the  signals  are  localized 
near  a  point  (0,yfl)  6  G  and  act  during  a  time  interval  T  <  t  <  0.  The 
resulting  wave  field  u(t,X)  is  then  characterized  by  its  initial  values 
u(0,X),  Dtu(0,X)  in  G.  To  make  explicit  their  dependence  on  y„  the 
initial  values  will  be  assumed  to  have  the  form 


(6.1) 


u(0,X)  »  f(X,y0)  =  f0(x,y-yfl)  , 
Dtu(0,X)  -  g(X,y0)  =  g0(x,y-y0) 


for  all  X  *  (x,y)  €  G  where  y 0  >  0, 

(6.2)  f0  €  Lj,COm(G)  ,  g0  €  Lj0m(G)  , 


and  f(X,yQ)  =  g(X,yQ)  =  0  for  (x,y-yQ)  G.  Note  that  for  y„  >  0  one  has 
f(’,y0)  €  D(Al/z),  g(*,y0)  €  L2 (G)  and  hence  u(t,X)  is  a  solution  wFE  in 
G.  The  functions  f(*,y0),  g(*»y„)  will  also  be  used  as  initial  values 
for  free  waves  in  R2  and  for  wave  fields  in  the  degenerate  grating 
domain  Rg .  In  each  case  the  domain  under  consideration  will  be  clear 


from  the  context  or  will  be  stated  explicitly.  For  brevity,  the 
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coordinate  y0  will  be  suppressed  except  in  places  where  the  y5 -dependence 
is  under  discussion. 

The  Signal  Wave  Field.  In  the  absence  of  a  diffraction  grating 
the  initial  state  f0,  g0  will  generate  a  signal  wave  field  us(t,X)  in  R2 . 
The  first  derivatives  of  ug(t,X)  have  asymptotic  wave  functions  (3, 
Theorem  2.10] 

(6.3)  Dk  ug(t,X)  -  r  1/2  sk(r-t,9)  +  o(l),  k  -  0,1,2  , 
where  the  waveforms  sk(T,0)  e  L2  (R  x  [-tt,tt])  are  given  by 

f  *00  x 

(6.4)  so(T,0)  -  Re  *iTU)  **#  (w  cos  sin  9)  (~iw) 1/2  dwj  , 

Sj(t,9)  -  -so(t,0)  cos  0,  s2(t,0)  *  -so(t,0)  a  in  0  and  the  terms  o(l)  -«•  0 
in  L2(R2)  when  t  -*■  «.  The  function 

(6.5)  fi0(P)  -  4>  g0(P)  -  iu(P)  <5  f0(P) 

where  $  denotes  the  Fourier  transform  in  L2(R2).  In  particular,  the 
Fourier  transform  in  L2  (R  *  [— tt. tt ] )  of  the  signal  waveform  so(t.0)  is 

(6.6)  % o  (to. 0)  *  j  (-iu>)  1/2  h0(u)  cos  0, u  sin  0)  . 

It  can  be  verified  that  if  f0  and  g0  are  real-valued  then  sq  (~  u>»  0) 

*  s0(o),9)  and  hence  (6.6)  generates  a  real-valued  signal. 

When  y0  is  large  the  signal  arriving  at  the  grating  surface  is 
described  by  the  signal  waveform  so(t,0)  through  (6.3).  The  problem  of 
signal  design  is  to  construct  a  source  or  "transmitter"  whose  waveform 
so(t,0)  approximates  a  prescribed  function.  The  solution  of  this  problem 
is  the  task  of  the  transmitter  design  engineer. 
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The  Echo  Wave  Fields.  In  Che  presence  of  a  diffraction  grating 
with  domain  G  the  initial  state  f,g  will  generate  a  total  wave  field 
u(t,X)  whose  asymptotic  behavior  for  t  -*• +«  is  described  by  Theorem  4.3. 
In  particular, 

(6.7)  D 0  u(t,X)  «  r~l/2  F0(r-t,6)  +  o(l)  in  L2(G),  t  ^ +®  , 

where  F0(T,8)  €  L2(R  *  [0, it] )  is  defined  by 


(6.8)  Fo(x,0)  =*  Re 


TTrfp72" 


.illD  c 


h  (id  cos  0,(d  sin  0)(-i(D)^2  doj} 


and 

(6.9)  h_(P)  =  g_(P)  -  iiD(P)  f_(P)  . 

The  echo  wave  field  ug(t,X)  is  defined  by 

(6.10)  ue(t,X)  -  u(t,X)  -  ug(t,X)  ,  t  >  0  ,  X  €  G  . 

Thus  the  echo  is  described  for  large  t  by 

(6.11)  D0ue(t,X)  -  r  1/2  e0  (r-t,  0)  +  o(l)  in  L2  (G) ,  t  -*•+<»  , 
where  e„  ■  F0  -  sQ  S  L2  (R  x  [ 0,  tt]  )  is  given  by 

(6.12)  eo(x,0)  *  Re  |fl  e^Ttii  h®c(w  cos  0,u)  sin  0)  (-ico)  dm 

with 

(6.13)  h®C(P)  -  h_(P)  -  h(P)  -  gfC(P)  -  iMP)  fSC(P)  . 
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The  last  functions  can  be  written  in  terms  of  the  R-B  diffracted  plane 
waves 


(6.14)  liL(X.P)  "  /nC(X,P)  +  ^fC(X,P) 

as  (5,  (2.24)] 


(6.15) 


g!C(P) 


l^SC  (X,P)  g(X)  dx 


Aap 

with  the  analogous  representation  for  f 

The  Echoes  of  Signals  from  Remote  Sources.  Equations  (6.12)  - 
(6.15)  provide  a  construction  of  the  echo  due  to  an  arbitrary  distribu¬ 
tion  of  sources.  The  principal  goal  of  this  section  is  to  determine  how 
this  construction  may  be  simplified  when  the  sources  are  far  from  the 
grating;  i.e.,  y„  ->•  <*>.  To  this  end  recall  the  decomposition  of  Lemma 
5.1.  Substituting  equation  (5.18)  in  (6.15)  gives 


m 


(6.16)  g_  (p,q)  “27 T  l  c^(p,q)  Kp^.-q^)  +  Pm>Q(p,q)  ,  (p,q)  <=  <?m  n  , 

&“-n  ’  ’ 


where  g  «  $  g  is  the  Fourier  transform  in  L2(R2)  and 


(6.17)  P(p,q)  =  P(p,q;g)  “  0  (X,p,q)  g(X)  dX  ,  (p,q)  €  R?  -  E  . 

'G 

Note  that  if  the  unitary  operator  R  :  L2(R2)  -►  L2(R2)  is  defined  by 

(6.18)  R  f(x,y)  -  f (x, -y) 


then 
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(6.19) 


(4>  R  f)(p,q)  =  (R  f )  (p,  q)  *  f(p,-q)  • 


Hence  (6.16)  implies  that  for  all  (p,q)  €  0^  Q  one  has 


(6.20) 


gfC(P) 


m 


2 "TT  l  c~(P)  (R  g)(P£)  +  Pm>n(P) 


A*-n 


-2*  i  c^cp)  « 

£=-m 


£  (S(R  §)_  ,  n  ,)  n(P)  +  Pm  n(P) 

n  ui  1  Jc  y  n-  Jo  in  y  n  ni  y  n 

Jo=-m 


-(SR  g)  _(P)  +  Pm  n(P) 
m,  n  m, n 


'  g£ 

by  Theorems  5.4  and  5.5.  Proceeding  in  the  same  way  with  w(P)  f_  (P) 
and  recalling  that  co(P)  *  u)(P^)  one  finds 


(6.21) 


hSC(P)  *  (SR  h)(P)  +  p(P;h)  ,  P  S  Rq  -  E  . 


where 


(6.22)  p(P;h)  -  P(P;g)  -  iw(P)  p(P;f)  . 

The  estimate  (5.19)  of  Lemma  5.1  clearly  implies  that 
p(P;h(*,y0))  -*•  0  when  y0  -*■  ",  uniformly  for  P  in  any  compact  subset  of 
Rq  -  E.  This  result  is  not  strong  enough  to  yield  a  corresponding 
estimate  of  the  echo  waveform  e0(x,8)  defined  by  (6.12)  and  it  is 
natural  to  conjecture  that  p(’;h(*,y0))  0  in  L2(R*)  when  y0  -*•  «. 

Unfortunately,  if  one  assumes  only  that  A(G)  admits  no  surface  waves 
then  this  property  does  not  follow  from  the  results  of  [5,6]  because  no 
information  was  obtained  concerning  the  behavior  of  i]/_(X,P)  for  P  near 
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the  exceptional  set  E.  However,  in  those  cases  where  the  analytic 
continuation  of  the  resolvent  of  Ap  has  no  singularities  on  o(Ap)  (i.e., 
Zp  n  a(Ap)  ■  <f>  for  every  p  6  (-1/2, 1/2])  the  limiting  absorption  theorem 
[5,  Corollary  4.17]  is  valid  on  all  of  a(Ap)  (see  [5,  Theorem  4.15])  and 
[5,  Theorem  6.1]  can  be  improved  to  state  that  t]>+(*,P)  exists  and 
P  -*■  iJj+(*,P)  S  L2,JloC(A,G)  is  continuous  for  all  P  €  R$".  This  improvement 
of  [5,  Theorem  6.1]  implies 

Theorem  6.1.  Let  A(G)  have  no  surface  waves  and,  in  addition, 
assume  that 

(6.23)  I  n  a(Ap)  =  <p  for  all  p  e  (-1/2, 1/2]  . 

Then  for  every  g0  6  L2  (G)  one  has 

(6.24)  §fC(*.y0)  “SR  g(,,y0)  +  o(l)  in  L2(R*),  yfl  -*•  °°  . 
Similarly,  for  all  f0  €  L*(G)  one  has 

(6.25)  w(*)  ffC(*,y0)  =  u)(0  S  R  f (•  ,y0)  +  o(l)  in  L2(R*),  yQ  ■+■  <»  . 

The  proof  of  Theorem  6.1  will  be  based  on  the  following  extension 
of  Lemma  5.1. 

Lemma  6.2.  Under  the  hypotheses  of  Theorem  6.1,  for  every 
compact  set  KC  s[  and  every  r'  >  r  >  h  there  is  a  constant  C  =  C(K,r,r') 
such  that 

(6.26)  |o+(X,P)|  <  C  for  all  X  e  R*,  and  P  e  K  . 

Proof  of  Lemma  6.2.  It  clearly  suffices  to  prove  the  lemma  for 
the  case  K  -  0^  q.  On  examining  the  proof  of  Lemma  5.1  one  finds  that 
the  continuity  of  P  +  ip+(*  ,P)  for  all  P  €  Rf  implies  that  (5.24)  holds 
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for  K  *  0  .  Moreover,  (5.25)  holds  for  all  P  €  0  with  u ,  «  u  »  0. 

m,n  in,  n  +  — 

Thus  (6.26)  follows  from  (5.27)  with  y(K)  ■  0. 

Proof  of  Theorem  6.1.  Note  first  that  if  the  translation 

operator  T  :  L2 (G)  ■*  L2 (G)  is  defined  for  each  y0  >  0  by  T  g. 

y0  “  Yo  0 

“  g(*,y0)  then  (6.24)  is  equivalent  to  the  statement  that 
(6.27)  s-lim  ($>_  -  i  -  S  R  $)T  *  0  . 

v  •+«>  y  o 

y  o 

Moreover,  the  family  of  operators  in  (6.27)  is  uniformly  bounded  for  all 

y0  >  0.  Hence  by  a  familiar  density  argument  (cf.  [3,  proof  of  Theorem 

2.6])  it  will  suffice  to  establish  (6.24)  for  all  g0  in  a  dense  subset 

00 

of  1.2 (G) .  The  set  C0(G)  will  be  chosen  for  this  purpose.  Thus  the 
proof  will  be  completed  by  showing  that  if  g0  €  C*(G)  and 


P(P,g(-,yo))  -  gf  (P,y„)  -  S  R  g( • ,y0) (p) 


(6.28) 


[  °_(x, 


P)  g(X,y0)  dX  ,  P  €  Rg  , 


then 


(6.29)  lim  [  |p(P,g(',yo)) | 2  dP  -  0  . 

y0-^Ro 

To  prove  (6.29)  it  will  be  convenient  to  decompose  R^  as  the  disjoint 
union 

(6.30)  Rq  -  D(y)  u  (D’(y)  n  E^)  u  (D’(y)  -  ^ ) 


where 
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(6.31) 


D(Y)  =  Rj  n  {P  :  |P|  >  y>  , 

D'(Y)  -  Ro  n  (P  :  [P |  <  Y>  ,  and 

E6  -  R2  n  {p  :  dist  (P,  )  <  6}  . 


With  this  notation  the  integral  in  (6.29)  can  be  written 


(6.32)  | p (P, g( * , y0) >  |  2  dP  =  Ii(Y.yo)  +  I2(Y»<5»y0)  +  MY.S.yo) 

JRo 


where 


(6.33)  Ii(Y»y0)  a  |p(P.g(*,y0) |2  dP  , 

JD(Y) 


(6.34)  I2(Y,'5,y0)  =  |p(P,g(*,y0)>  I  2  dP  »  and 

jD'(Y  )HE5 


(6.35)  I3(Y.<S,y  )  *  I p(P,g(- ,y0))  I  2  dP  . 

V(Y)-e5 


To  estimate  I1(Y,yo)  note  that  (5.18)  implies  that  (A  +  |p|2)  a_(X,P)  =  0 
for  all  X  €  G,  P  €  R2.  Thus  integrating  by  parts  in  (6.28)  gives 


p(p,g(*,y0))  »  -j p|  Act  (X,P)  g(X,y0)  dX 

JG 

(6.36) 

-  -|Pj"2  a  (X,P)  Ag(X,y0)  dX 
G 

-  -|P|2  (Ag(-,y0)2,SC(P)  -  S  R(Ag(-,y0))"(p)} 


-jP|"2  {(Ag(-,y0)f(P)  -  (Ag(*  ,y0))A  (P)  -  S  R(Ag(  •  ,y0)  f(P)  } 
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On  squaring  (6.36),  integrating  over  D(y)  and  using  the  inequality 
\zl  +  z2  +  z3j2  <  4  ( I  z  j.  | 2  +  |  z  2 1  2  +  |  z  3 1 2  >  one  finds 

Ii(Y,y0>  5  [  !pr4  |CAg>2  -  (Ag)~  -  S  R(Ag)  ^  |  2  dP 
;D(y) 

(6.37) 

<  4y~4  ( 1  <l>_  ( Ag  (*,y0))«2  +  »*(Ag(-,y0)l2  +  IS  R(Ag(-,y0))l 2) 

<  12y'“  8Ag(-,y0)l2  -  12y"“  IA8oIl2(G) 

for  all  y0  >  0.  In  particular,  Ij(Y,y0)  is  small  for  large  y, 
uniformly  in  y0  >  0. 

Now  consider  I2(Yf<$>y<)).  Lemma  6.2  and  equation  (6.28)  imply 
that  for  all  P  6  D '  (y)  one  has 


(6.38) 


|p(P,g(*.y0) 1  <  C(D'(Y),r,r')  J  |g0(x,y-yo)|  dx  dy 

-  Cj(Y,r,r')  [  |g,(X)|  dX  »  C2 (g0 .y.r, r ')  . 

*  c. 


Combining  this  and  (6.34)  gives 

(6.39)  I2(Y,6,y0)  <  C2(g0,Y,r,r’)  |D'(Y)  n  )  , 

for  all  y0  >  0,  where  |m|  denotes  the  Lebesgue  measure  of  a  set  M  C  R2 . 
Finally,  note  that  Lemma  5.1  implies  that  p(P,g(*,y0))  -►  0  when 
y„  -*•  «,  uniformly  for  P  6  D'(y)  -  Ej,  when  y  >  0  and  6  >  0  are  fixed. 
Thus 


(6.40)  lim  I3(Y>6,y0)  *  0  ,  y  and  6  fixed  . 

v  -K® 
y  0 
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To  complete  the  proof  of  (6.29)  let  e  >  0  be  given  and  use  (6.37) 
to  choose  a  y  =  y0  *  Yo(£«g0)  >  0  such  that  Ii(y,yo)  <  e/3.  Next  use 
(6.39)  with  Y  ■  Y0(£*ga)  fixed  to  choose  6  ■  60  *  6o(e.g0)  >  0  so  small 
that  l2(Yo»5o»y„)  <  e/3.  Both  of  these  estimates  hold  uniformly  for  all 
y0  >  0  .  Finally,  choose  Y0  -  Y0(e,g0)  so  large  that  I3  (Y0  ><5,,  ,y„)  <  e/3 
for  all  y0  >  Y0.  This  is  possible  by  (6.40).  With  these  choices  (6.32) 
implies  that 


(6.41) 


|p(P.g(*,y0))|2  dP  <  e  for  all  y.  > 


Ve,  80) 


which  proves  (6.29)  and  therefore  (6.24).  Finally,  to  prove  (6.25)  one 
notes  that  if  f(*,y0)  6  L*  (G)  then  ui(P)  f_(P,y0)  S  L2(R*)  and  the 
preceding  argument  can  be  applied  to  this  function.  This  completes  the 
proof  of  Theorem  6.1. 

An  Estimate  of  the  Echo  Waveform.  Under  the  hypotheses  of 
Theorem  6.1  one  has  the  estimate 

hfc(*,y0)  -  gfc(*.y0>  -  iw(0  ffc(-,y0) 

(6.42) 

•SR  g(*,y„)  -  ioj(*)  S  R  f(*,y0)  +  0(1) 

•SR  h(‘,y0)  +  o(l) 

where  o(l)  0  in  L2(R*)  when  yfl  -►  ®.  Moreover,  the  mapping  hfc  €  L2(R*) 

■*  e0  €  L2  (R  x  [ 0, tt ] )  defined  by  (6.12)  is  bounded  with  bound  1  [3,  (2.84)]. 
It  follows  that 


.  00 

V/W1T17S  I  eiTu)(S  R  h)  (wcos  9,u)  sin  0,yo)  (-iw)  ^  dw|+ o(l) 


(6.43)  e0 (x, 9)  ■  Re 


o 


where  o(l)  ■»  0  in  L2(R  x  [O.tt])  when  y0  -*■  ».  Now 


(6.44) 


h(p.q»y0)  ■  §(p»q»yo>  -  iw(p.q)  f(p»q.y0) 

.  eiqy°[g0(p,q)  -  iu(p,q)  f0(p,q)] 
-  h,(p,,) 


and  hence  by  (6.6) 

(-iu>)  R  h(ucos  6,0)sin  0,yo) 
(6.45) 


(-id>)  ^2  h(u)  cos  0 ,  —to  sin  9,y0) 
e~ia)y°  sin  ®(-id))  l^2  h„  (to  cos  0,-w  sin  9) 
2  8<fol,_e)  . 


Combining  (6.43)  and  (6.45)  gives 


(6.46) 


eo(T,0)  -  Re 


2 1  1/2  |  eixw  e-iwyo  sin  9 


S(u>)  sn  (u),-0)du)f  +  o(l) 


Thus  under  the  hypotheses  of  Theorem  6.1  the  echo  waveform  is  determined 
by  the  signal  waveform,  the  S-matrix  for  the  grating  and  the  range 
parameter  y0,  with  an  error  that  tends  to  zero  in  energy  when  y0  -►  °°. 

Pulsed  Beam  Signals.  For  many  applications  it  is  desirable  to 
have  a  transmitter  whose  waveform  so(t,0)  is  sharply  limited  in  both 
direction  and  frequency.  The  relation  (6.6)  shows  that  this  could  be 
achieved  by  choosing  f.  and  g.  such  that  supp  h  «  RK  where  K  c  0  and 
m  and  n  are  suitably  chosen.  Of  course,  this  condition  cannot  be  satis¬ 
fied  with  sources  that  are  confined  to  a  compact  set,  since  hQ (P)  is  then 
analytic.  However,  it  may  be  possible  to  choose  f0  and  g0  such  that 


» 
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(6.47) 


h.(p,q)  -  a(p,q)  +  b(p,q) 


where 


(6.48) 

(6.49) 


supp  a  *  RK  C  RO  , 
in,  n 


s3(to,9)  ■  ^(-iw)  a  (to  cos  9,uisin  9) 


defines  the  desired  waveform  s“  and 


(6.50) 


lb,L2(R2)  <  £  ' 


If  this  transmitter  design  problem  has  been  solved  then  the  corresponding 
echoes  will  satisfy 


(6.51) 


e0 (t,9)  -  ef(x,9)  +  Oj  +  o 


where 


(6.52)  ej(t,9)  -  Re  j  | 


1/2 


eiT0)  e-iioyo  Sin  9  s  A  sa 


S (to)  8*(u,-9)duJ 


while 


(6.53) 


i  o.l 


i'l2(Rx[0,tt])  <  £  for  311  >  0  , 


and 


(6‘54)  ,O*SL2(BX[0.ir])  "°  ’ 

70 

Angular  Dispersion  of  Echoes  from  Gratings.  The  notation 


(6.55) 


r.  -  {P  -  (oicos  9,u  sin  9)  :  go  >  0  and  oto  <  9  <  6«) 


will  be  used  to  denote  the  smallest  sector  such  that  RjC  -n  <  l  <  m. 
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The  hypothesis  K  »  Ka  C  0^  ^  implies  that  the  sectors  are  disjoint  and 

(6.56)  Q  r£  C  R*  . 

£■- n 

Moreover,  (6.48),  (6.49)  and  (6.52)  and  Theorem  5.4  imply  that  one  has 

m 

(6.57)  supp  r~l/z  e^(r-t,0)  C  y 

£*>-n 

for  all  t  >  0.  Thus,  apart  from  the  error  terms  in  (6.51),  the  echo 
waveform  is  concentrated  in  the  sectors  I^.  Note  that  in  the  case  of  a 
degenerate  grating  with  Neumann  (resp.,  Dirichlet)  boundary  condition 
one  has  S  *  1  (resp.,  S  ■  -1)  and  hence  r  e®(r-t,0)  «  ±r~1/,z  s^(r-t,0) 
has  support  in  rQ.  This  is  a  well-known  property  of  the  specular 
reflection  of  a  beam  by  a  plane.  In  the  case  of  a  non-degenerate  grating, 

/v 

where  S  /  ±1,  one  has  only  (6.57)  and  secondary  reflected  beams  will 
appear  in  the  sectors  I^,  i  i  0.  Their  waveforms  can  be  calculated 
explicitly  using  (6.52)  and  (5.87).  They  are  distortions  of  the  signal 
waveform  so(t,0)  whose  forms  are  determined  by  the  scattering  coeffi¬ 
cients  c^(oj  cos  0,oj  sin  0) .  This  phenomenon  of  the  angular  dispersion  of 
pulsed  beams  by  diffraction  gratings  is  the  counterpart  for  transient 
wavefields  of  the  phenomenon  of  the  diffraction  of  monochromatic  beams 
into  the  higher  order  grating  directions. 
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